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Role of structural fluctuations and environmental 
noise in the eiectron/hole separation kinetics at or¬ 
ganic polymer bulk-heterojunction interfaces. 

Eric R. Bittner and Allen Kelley ^ 


We investigate the electronic dynamics of model organic photovoltaic (OPV) system consisting of 
polyphenylene vinylene (PPV) oligomers and [6,6]-phenyl C61-butyric acid methylester (PCBM) 
blend using a mixed molecular mechanics/quantum mechanics (MM/QM) approach. Using a 
heuristic model that connects energy gap fluctuations to the average electronic couplings and 
decoherence times, we provide an estimate of the state-to-state internal conversion rates within 
the manifold of the lowest few electronic excitations. We find that the lowest few excited states of a 
model interface are rapidly mixed by C=C bond fluctuations such that the system can sample both 
intermolecular charge-transfer and charge-separated electronic configurations on a time scale of 
20fs. Our simulations support an emerging picture of carrier generation in OPV systems in which 
interfacial electronic states can rapidly decay into charge-separated and current producing states 
via coupling to vibronic degrees of freedom. 


1 Introduction. 

The power conversion efficiencies of highly optimized organic 
polymer-based photovoltaic (OPV) cells exceeds 10% under stan¬ 
dard solar illuminationlil with reports of efficiencies as high 
as 12% in multi-junction OPVs. This boost in efficiency indi¬ 
cates that mobile charge carriers can be generated efficiently in 
well-optimized organic heterostructures; however, the underly¬ 
ing mechanism for converting highly-bound Frenkel excitons into 
mobile and asymptotically free photocarriers remains elusive in 
spite of vigorous, multidisciplinary research activity.HEl Much of 
the difficulty in developing a comprehensive picture stems from 
our lack of detailed understanding of the electronic energy states 
at the interface between donor and acceptor materials and how 
these states are influenced by small-but significant-fluctuations 
in the molecular structure. 

The generic photophysical pathways that underlie the genera¬ 
tion of mobile charge carriers in an OPV solar cell are sketched 
in Fig. The absorption of a photon by the material produces 
a TT — TT* excitation (exciton) within the bulk (a) that can mi¬ 
grate and diffuse via Forster energy transfer processes. Once 
in close proximity to a bulk-heterojunction interface, energetic 
off-sets between the respective HOMO and LUMO levels of ad- 
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Fig. 1 Jablonski diagram showing relative energetic positions of excited 
states in a bulk-heterojunction system. The discrete states on the left 
are assumed to be localised excitons or charge-transfer states while the 
polaron states are assumed to be a quasi-continuum of mobile 
charge-separated states within the bulk. The latter of these are 
presumed to be responsible for producing photocurrent from a given 
device. 


jacent donor and acceptor molecules provide the necessary driv¬ 
ing force to separate an exciton into a localized charge transfer 
(CT) state (b) which typically lies 0.25 to 0.4 eV lower in en¬ 
ergy. Alternatively, an exciton may dissociate directly via tunnel¬ 
ing into charge-separated (CS) or polaron states (c) which may 
subsequently evolve to contribute to the photocurrent or undergo 
geminate or non-geminate recombination to form CT states (d). 
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Fig. 2 Snapshot of the MD simulation cell containing 50 PCMB 
molecules and 25 PPV oligomers following equilibration at 100 K and 
Ibar of pressure. The red and blue highlighted molecules denote the 
TT-active units in our system. 


We distinguish CT states from CS states by whether or not the 
donor and acceptor species are in direct contact (CT) or sepa¬ 
rated by one or more intermediate molecules (CS). 

Ultrafast spectroscopic measurements on OPV systems have re¬ 
ported that charged photoexcitations are generated on < 100-fs 
timescales |2|i2 | i5 | i9j [22l however, full charge separation to pro¬ 
duce free photocarriers is expected to be energetically expen¬ 
sive given the strong Coulombic attraction between electrons and 
holes due to the low dielectric constant in molecular semiconduc¬ 
tors. Nonetheless, experiments by Gelinas et ah, in which Stark- 
effect signatures in transient absorption spectra were analysed to 
probe the local electric field as charge separation proceeds, indi¬ 
cate that electrons and holes separate by as much as 40A over the 
first 100 fs and evolve further on picosecond timescales to pro¬ 
duce unbound and hence freely mobile charge pairsEI. Concur¬ 
rently, transient resonance-Raman measurements by Provencher 
et al demonstrate clear polaronic vibrational signatures on sub- 
100-fs on the polymer backbone, with very limited molecular 
reorganization or vibrational relaxation following the ultrafast 
stepl^. Such rapid through-space charge transfer between exci- 
tons on the polymer backbone and acceptors across the hetero¬ 
junction would be difficult to rationalize within Marcus theory 
using a localised basis without invoking unphysical distance de¬ 
pendence of tunnelling rate constantsES and appear to be a com¬ 
mon feature of organic polymer bulk heterojunction systems. 

Polymer microstructural probes have revealed general relation¬ 
ships between disorder, aggregation and electronic properties in 
polymeric semiconductorsEH. Moreover, aggregation (ordering) 
can be perturbed by varying the blend-ratio and composition of 
donor and acceptor polymersEl. On one hand, energetic dis¬ 
order at the interface would provide a free-energy gradient for 
localised charge-transfer states to escape to the asymptotic re¬ 
gions. In essence, the localised polarons in the interfacial region 
could escape into bands of highly mobile polarons away from 
the heterojunction regionES On the other hand, energetic dis¬ 


order in the regions away from the interface would provide an 
entropic driving force by increasing the density of localised po- 
laron states away from the interfacial region, allowing the po¬ 
larons to hop or diffuse away from the interface before recombi¬ 
nation could take placeESI. it has also been suggested that in poly- 
mer/fullerene blends, interfacial exciton fission is facilitated by 
charge-delocalisation along the interface which provide a lower 
barrier for fission with the excess energy provided by thermally- 
hot vibronic dynamicsEZl. Finally, a report by Bakulin et al in¬ 
dicates that when relaxed charge-transfer-excitons are pushed 
with an infrared pulse, they increase photocurrent via delocalised 
states rather than by energy gradient-driven hoppingEZl. 

Recent MD simulations by Fu et al of a model 
Squarene/Fullerene OPV cell indicate that the degree of disor¬ 
der at the interface directly affects couplings and hence golden- 
rule transition rates between the ground and excited states. The 
disorder in the system is at least partly introduced by thermal an¬ 
nealing and the interactions between the Squaraine and Fullerene 
layers. Their simulations indicate that even at BOOK, the ther¬ 
mal motions of the molecules at the interface can be quite pro¬ 
found and the degree of disorder inherent around the interface 
can greatly affect the formation of electron/hole pairs .ESI 

Finally, we recently presented a fully quantum dynamical 
model of photo-induced charge-fission at a polymeric type-II het¬ 
erojunction interface .ESI Our model supposes that the energy level 
fluctuations due to bulk atomic motions create resonant condi¬ 
tions for coherent separation of electrons and holes via long- 
range tunnelling. Simulations based upon lattice models reveal 
that such fluctuations lead to strong quantum mechanical cou¬ 
pling between excitonic states produced near the interface and 
unbound electron/hole scattering states resulting in a strong en¬ 
hancement of the decay rate of photo excitations into unbound 
polaronic states. 

In this work, we employ the our model in conjunction with an 
atomistic quantum/molecular dynamics study of a model donor- 
acceptor blend to provide an estimate of the interstate transition 
rates with the excited-state manifold. Our simulations combine a 
molecular dynamics description of the atomic motions coupled to 
an on the fly description of the excited state 7r-electronic structure. 
By analysing the energy gap fluctuations between excited states, 
we provide a robust estimate of both the interstate electronic cou¬ 
plings, decoherence times, and transition rates. We begin with a 
brief overview of our model and then describe our results. 

1.1 Heuristic model 

A simple model for considering the role of the environment can be 
developed as followsESISni. Consider a two-state quantum system 
with coupling A in which the energy gap A{t) fluctuates in time 
about its average and (A(?)A(0)) = In a two-state basis of 
{|0), |1)} the Hamiltonian for this can be written as 

H=^a^ + ?iax ( 1 ) 
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where 6)^ are Pauli matrices. Note that Eq. [^can be transformed 
such that fluctuations are in the off-diagonal coupling 


H — —dz-\- {t)dx 


( 2 ) 


where = A +A and 3V(t) = 0. The fluctuations in the electronic 
energy levels are attributed to thermal and bond-vibrational mo¬ 
tions of polymer chains and can be related to the spectral density 
S{co) via 


V^ = 5V2{t) = jy^S{(0). (3) 

Averaging over the environmental noise, we can write the average 
energy gap as hCl = y^AjTW with eigenstates 

|i/r+) = cos0|O)+ sin0|l) (4) 

|i/r_) = — sin0|O)+ cos0|l) (5) 

where tan20 = |t^|/Ao defines mixing angle between original kets. 
Consequently, by analysing energy gap fluctuations, we can ob¬ 
tain an estimate of both the coupling between states as well as 
transition rates. 


1.2 Noise-averaged kinetics 


To estimate the average transition rate between states, we first 
write the equations of motion for the reduced density matrix for 
a two level system coupled to a dissipative environment, 

/ - , X 1 

Pll = - tV(P 21 - Pi 2 ) - — Pn 

n T\ 

i - ^ X 1 

P22 = TzV[P2l - Pl2) -P22 

in T 2 

Pl2 = Tt^(P22-Pll) - 7^Pl2 + ^Pl2 (6) 

n n 

Pll = ^^(P22-Pll) - ^P21-^^21 

where we have introduced ti and T 2 as the natural lifetimes of 
each state and as the decoherence time for the quantum su¬ 
perposition, which in turn, can be related to the spectral density 
via =V/fi. Taking to be short compared to the lifetimes of 

each state, we can write the population of the initial state as 


pii(0 = exp 



(7) 




(c) 


Fig. 3 (a)Single Cl (SCI) energies for the first 50 fs following excitation 
to the lowest SCI state, (b) Histogram distribution of the 5 lowest 
excitation energy levels over a 50ps simulation. Throughout the 
simulation, the lowest lying state (in blue) remained populated and 
varied in character from excitonic to charge-separated. Labels A-E refer 
to the configurations shown in Fig. 5. The color of each component 
corresponds to specific SCI eigenvalue depicted in (a), (c) Histogram of 
energy gaps between the first excited and higher-lying excited states. 
Coldi 1 —y 2, Blue! 1 —y 3, Qreeni 1 —y 4, Red! 1 —y 5. 


where k is average state-to-state transition rate. Integrating this 
over time, we obtain an equation of the form 



suggesting a form for the exact solution of Eqs[^ Taking the 
Laplace transform of Eqs|^ and assuming that our initial popu¬ 
lation is in state 1 (pii(O) = 1), Eqsj^become a series of algebraic 
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equations 


-1 = - W{p 2 i -Pn) - — Pii 

n Ti 

0 = TViPii-Pn) -P 22 

n T 2 

0 = 7 t/(p 22 -Pll) - —P 12 +^Pl 2 (9) 

h n 

r. ^ ~ ~ \ 1 ~ ~ 

0 = -7^(P22-Pll) - ;^P21-^P21 

n n 

which after a bit of algebra gives a rate constant of the form 


k 


Td 

{Td^olhY + \ 


( 10 ) 


The average rate vanishes in the limit of rapid decoherence {T^ 

0). This is the quantum Zeno effect whereby rapid quantum mea¬ 
surements on the system by the environment at a rate given by 
collapses the quantum superposition formed due to the in¬ 
teractions and thereby suppresses transitions between states. On 
the other hand, transient coherences can facilitate quantum tran¬ 
sitions between otherwise weakly coupled states. 

For the sake of connecting this to the photophysical dynamics 
of a OPV heterojunction, let us assume that one state corresponds 
to a CT state and the other corresponds to a CS state. When the 
fluctuations are weak, V <c A^, the original kets |Cr) and |CA) pro¬ 
vide a good zeroth description of actual eigenstates of the system 
and the coupling can be treated as a weak perturbation. How¬ 
ever, when the off-diagonal couplings are comparable to the av¬ 
erage gap, even low-lying CT states may be brought into strong 
coupling with the CS states. 

The implication of this heuristic model is that one can obtain 
the input to the rate equation (Eq. [T^ from mixed quantum 
classical simulations that take into account the excited state pop¬ 
ulations. Here we report on such simulations of a model OPV 
system consisting of a blend of fullerene and polyphenylene viny- 
lene oligomers as depicted in Fig. Our simulation method em¬ 
ploys an atomistic description of the nuclear dynamics described 
by a force-field that responds to changes in the local n electronic 
structure of a sub-set of molecules with in the simulation cell. We 
restrict the excited state population to the lowest tt — tt* excitation 
as to simulate the long-time fate of a singlet CT state prepared 
via either photoexcitation or charge recombination. By analysing 
the energy gaps between electronic adjacent states and the char¬ 
acter of the excited states in terms of electron-hole configuations 
we can deduce how small vibronic motions of the polymer chains 
modulate the electronic coupling and induce charge-separation. 


2 Methods. 


Our simulations employ a modified version of the TINKER molec¬ 
ular dynamics (MD) packageED in which the MMS^ZI intramolec¬ 
ular bonding parameters were allowed to vary with the local 
TT-electronic density as described by a Parisier-Parr-Pople (PPP) 
semi-empirical Hamiltonian .1^331 Specifically, we assume that the 
internal bond force constants, bond-lengths, bond angle, bend¬ 


ing potentials, and bond torsion parameters are linear functions 
of the local bond-order. We specifically chose 1 PCBM and 3 
nearby PPV oligomers to represent a model bulk-heterojunction 
in order to study the penetration of extended intramolecular elec¬ 
tronic states into the bulk region. The remaining molecules in 
the simulation were treated using purely classical force-field. At 
each step of the simulation, we compute the Hartree-Fock (HF) 
ground-state for the n system and use configuration interaction 
(singles) (SCI) to describe the lowest few n ^ 7f excitations. In- 
termolecular interactions within the active space were introduced 
via non-bonding Coulombic coupling terms and static dispersion 
interactions contained within the MM3 forcefield. 

Fig.[^shows a snapshot of the periodic simulation cell contain¬ 
ing 50 PCMB molecules and 25 PPV oligomers following equili¬ 
bration at lOOK and Ibar of pressure using classical molecular 
dynamics (MD) within the NPT ensembl^ The molecules sur¬ 
rounding the 4 TT-active molecules serve as a thermal bath and 
electronic excitations are confined to the 7r-active orbitals. In to¬ 
tal, our TT-active space included a total of 172 carbon Ip^ orbitals 
and we used a total of 10 occupied and 10 unoccupied orbitals 
to construct electron/hole configurations for the Cl calculations. 
During the equilibration steps, we assume the system to be in 
its electronic ground state, after which we excite the system to 
the first SCI excited state and allow the system to respond to 
the change in the electronic density within the adiabatic/Born- 
Oppenheimer approximation. It is important to note that the ex¬ 
cited state we prepare is not the state which carries the most 
oscillator strength to the ground nor do we account for non- 
adiabatic surface hopping-type transitions in our approach 
Our dynamics and simulations reflect the longer-time fate of the 
lowest-lying excited state populations as depicted in step (d) in 
Fig. 0 The combination of a classical MD forcefield with a semi- 
empirical description of a select few molecules within the system 
seems to be a suitable compromise between a fully ab initio ap¬ 
proach which would be limited to only a few molecules and short 
simulation times and a fully classical MD description which would 
neglect any transient changes in the local electronic density!^. In 
spite of the relative simplicity of our model, the simulations re¬ 
main quite formidable. 

3 Results 

Over the course of a 50 ps simulation, the lowest lying SCI exci¬ 
tation samples a variety of electronic configurations ranging from 
localised PCBM excitons to charge-separated and charge-transfer 
states with varying degrees of charge separation. In Fig. [^, we 
show the lowest few SCI excitation energies following excitation 
at ? = 0 fs to the lowest SCI state for one representative simu¬ 
lation. The labels refer to snap-shots taken at 50fs intervals to 
visualize the various electronic configurations sampled by our ap¬ 
proach in Fig. First, we note that following promotion to the 
lowest lying SCI at ? = 0 there is very little energetic reorganiza¬ 
tion or relaxation compared to the the overall thermal fluctua¬ 
tions that modulate the SCI eigenvalues. This can be rationalized 


*PCBM:Phenyl-C61-butyric acid methyl ester, PPV: Polyphenylene vinylene 
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(a) 0 fs (b) 50 fs (c) 100 fs 


0.0% 2.5% 


29 . 5 % 67 . 7 % 


61 . 2 % 29 . 8 % 





35 . 0 % 38 . 9 % 





(d) 150 fs (e) 100 fs (f) 150 fs 

Fig. 4 Sequence of states following initial excitation. In this particular simulation, the initial excited state population was localised on the PCBM. This 
population rapidly dissociates into charge-transfer and then charge-separated configurations within 100-150 fs. This scenario appears to be a 
ubiquitous feature of our simulations indicating that fluctuations alone can carry the system from a localized CT state to charge-separated states and 
eventually current-producing polaron states. 


since the excitation to the ;r* orbitals at any given time is largely 
delocalised over one or more molecular units and hence the av¬ 
erage electron-phonon coupling per C=C bond is quite small. 
The number of avoided crossings that occur between lowest ly¬ 
ing states is highly striking, signalling that the internal molecular 
dynamics even at lOOK is sufficient to bring these states into re¬ 
gions of strong electronic coupling. 

In Fig. we show a histogram of the energies of the lowest 5 
SCI energies accumulated over 50 ps of simulation time follow¬ 
ing promotion to the lowest SCI state. Focusing upon the lowest 
lying state (in blue), fluctuations due to molecular dynamics can 
account for nearly 0.25eV of inhomogeneous broadening of this 
state bringing it (and similarly for the other states) into regions of 
strong electronic coupling with nearby states which dramatically 
changes the overall electronic character of the state from purely 
charge-transfer to purely excitonic on a rapid time-scale. 

We next consider the origins of the energy fluctuations evi¬ 
denced in Fig.[^. While we only show a 200fs segment of a much 
longer 5Ops simulation, over this period one can see that the Cl 
energies are modulated with a time-period of about 20 fs. Indeed, 
the Fourier transform of the average Cl excitation energy reveals 
a series of peaks between 1400-1700 cm“^ which corresponds 
to the C=C stretching modes present in the polymer chains and 
fullerene molecules. We conclude that small-scale vibronic fluc¬ 
tuations in the molecular structures and orientations produce sig¬ 
nificant energetic overlap between different electronic configu¬ 
rations in organic polymer-fullerene heterojunction systems and 


speculate that this may be the origins of efficient charge separa¬ 
tion in such systems. 

A variety of states are produced by our model and we cate¬ 
gorize them as excitonic (electron/hole on the same molecular 
unit) occupying 50% of the states, charge-transfer (electron on 
PCMB/hole on PPV#1 (nearest to PCBM)) occupying 15% of 
the states, or charge separated (electron on PCBM/hole on PPV 
#2 or #3) occupying 25% of the states. These are shown in 
Fig. HI where we indicate the local charge-density of each C2pz 
orbital by a blue (positive) or red (negative) sphere scaled in pro¬ 
portion to the local charge on each site. 

In Fig. we show a sequences of snap shots taken every 50 
fs following the initial excitation. These figures reveal the highly 
dynamical nature of the lowest lying excitation in which the sys¬ 
tem samples excitonic configurations (a), charge-transfer config¬ 
urations (b,c), charge-separated states in which the electron and 
hole are separated by at least one polymer chain (d,e), and highly 
mixed charge-tranfer/excitonic states as in (e,f). The scenario 
depicted here is not an exceptional case within our simulation. 
Here, the initial excited state is clearly localized on the PCBM unit 
with a high degree of excitonic character. There is clearly some 
charge separation within the PCBM; however, both the electron 
and hole clearly reside on the PCBM. After 50 fs (c.f. Fig. 
structural fluctuations in part induced by changes in electronic 
density and in part by the thermal fluctuations of the surrounding 
medium bring this state nearly into resonance with the second SCI 
state inducing charge separation between the PCBM and a nearby 
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(d) 54 fs (e) 55 fs (0 56 fs 

Fig. 5 A 6 fs time history from 51 to 56 fs of the first excited Cl state is shown in Fig. This series resides around an avoided crossing between the 
first and second excited state. Showing how a CT state at the interface can evolve into exciton states farther into the bulk around one avoided crossing 
points. After the excited states diverge, the exciton in the bulk collapses back into an exciton localised on the PCBM. 


PPV oligomer. Further fluctuations bring this state into resonance 
with other charge separated states as shown in Fig. Ef. d, and 
f. Because our quantum subspace is restricted to the molecules 
shown in Fig|^ charge separation to more distant regions can not 
occur. 

Fig.|^ shows a histogram of the energy gaps between the first 
and higher-lying excited excited states of our system. According 
to our heuristic model, the mean (ju) and variance (cr^) of the gap 
distributions can be used as input to estimate the state-to-state 
rate for this system under the mapping that jU = and cr^ = dV^ 
as per Eq.[^ Furthermore, we take ^V/h as an estimate of 
the decoherence time. The results are shown in Table [T] The es¬ 
timated transition rates are consistent with the observation that 
the system rapidly samples a wide number possible configurations 
over the course of the molecular dynamics simulation. On aver¬ 
age, the state-to-state coupling of 70 meV is comparable to the 
average energy gaps between the lowest states. This brings the 
system into the regime of strong electronic coupling. However, 
larger couplings also imply shorter electronic decoherence times- 
which dramatically limit the the ability of charges to separate by 
tunnelling. In this case, any superposition Cl states resulting from 
the quantum mechanical time-evolution of the system would col¬ 
lapse to single eigenstate on a sub-10 fs time-scale due to the 
interaction with the vibronic degrees of freedom. Based upon our 
model and simulations, we estimate that within 20 fs, population 
in any low-lying electronic state is effectively mixed with nearby 
states simply due to the underlying vibrations of the lattice. 

Fig. illustrates how the electronic nature of the lowest lying 
6 I Journal Name, [year], [vol.],i-j^ 


excitation changes over a short time period. Shown here are the 
excess charge-densities at each atomic site taken at 1 fs intervals. 
At each 1 fs time-step, the electronic character changes wildly due 
to strong coupling between the electronic and vibrational degrees 
of freedom. 

Table 1 Estimated interstate transition rates from Eq.[j^and vibronic 
couplings. 


Transition 

Ao (eV) 

V (eV) 

Td (fs) 

(fs) 

1^2 

0.11 

0.06 

11.41 

21.8 

1^3 

0.18 

0.08 

8.66 

16.6 

1^4 

0.26 

0.08 

8.70 

16.7 


4 Conclusions 

We present here the results of hybrid quantum/classical simula¬ 
tions of the excited states of a model polymer/fullerene hetero¬ 
junction interface. Our results indicate that dynamical fluctua¬ 
tions due to both the response of the system to the initial excita¬ 
tion and to thermal noise can rapidly change the character of the 
lowest lying excited states from purely excitonic to purely charge 
separated over a time scale on the order of 100 fs. In many cases, 
an exciton localised on the fullerene will dissociate into a charge- 
transfer state with the hole (or electron) delocalised over multiple 
polymer units before localizing to form a charge-separated state. 
While the finite size of our system prevents further dissociation of 
the charges, the results are clearly suggestive that such interstate 
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crossing events driven by bond-fluctuations can efficiently sepa¬ 
rate charges to a distance to where their Coulombic attraction is 
comparable to the thermal energy. 

The nuclear motions most strongly coupled to the electronic de¬ 
grees of freedom correspond to C=C bond stretching modes im¬ 
plying that small changes in the local lengths have a dramatic role 
in modulating the electronic couplings between excited states. 
Since kT <c hco for these modes, they should be treated quan¬ 
tum mechanically rather than as classical motions. Generally 
speaking, including quantum zero-point effects into calculation 
of Fermi golden-rule rates leads to slower transition rates than 
those computed using classical correlation functions which im¬ 
plies that the values estimated here are the upper bounds for 
the actual transition rates. The results presented here corrobo¬ 
rate recent ultrafast experimental evidence suggesting that free 
polarons can form on ultrafast timescales (sub lOOfs) and un¬ 
derscore the dynamical nature of the bulk-heterojunction inter- 

face. l4l718il7|22J24|29|^39| 
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